Simulating Dense QCD Matter with Ultracold Atomic Boson-Fermion Mixtures 
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We delineate, as an analog of two flavor dense quark matter, the phase structure of a many- 
body mixture of atomic bosons and fermions in two internal states with a tunable boson-fermion 
attraction. The bosons b correspond to diquarks, and the fermions / to unpaired quarks. For 
weak b-f attraction the system is a mixture of a Bose- Einstein condensate and degenerate fermions, 
while for strong attraction composite b-f fermions N, analogs of the nucleon, are formed, which 
are superfluid due to the A'^-A'^ attraction in the spin-singlet channel. We determine the symmetry 
breaking patterns at finite temperature as a function of the b-f coupling strength, and relate the 
phase diagram to that of dense QCD. 
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Ultracold atomic systems and high density QCD mat- 
ter, although differing by some 20 orders of magnitude in 
energy scales, share certain analogous physical aspects, 
e.g., crossovers from Bose-Einstein condensation (EEC) 
to BCS 0, . Motivated by phenomenological studies of 
QCD that indicate a strong spin-singlet diquark correla- 
tion inside the nucleon [3] , we focus here on modeling the 
transition from color superconducting two-flavor quark 
matter (2SC) at high density to superfluid hadronic mat- 
ter at low density in terms of a boson-fermion system, in 
which small size diquarks are the bosons, unpaired quarks 
the fermions, and the extended nucleons are regarded as 
composite boson-fermion particles. Recent advances in 
atomic physics have made it possible indeed to realize 
such systems in the laboratory. In particular, tuning the 
atomic interaction via a Feshbach resonance allows for- 
mation of heteronuclear molecules, as recently observed 
in a mixture of *^Rb and '^''K atomic vapors in a 3D 
optical lattice [4| , and in an optical dipole trap [5] . 

The analogy is incomplete however. The gluonic at- 
traction in QCD is a function of the baryon density, and 
thus tuning the coupling strength at fixed density is not 
possible in dense matter; furthermore, chiral symmetry 
breaking plays an important role in the QCD transition 
0. With these reservations in mind, we suggest that 
fuller understanding, both theoretical and experimental, 
of the boson-fermion mixture, as well as a mixture of 
three species of atomic fermions, as discussed in Q, can 
reveal properties of high density QCD not readily observ- 
able in laboratory experiments. 

We first delineate possible phase structures that can 
be realized in a mixture of single-component bosons (&) 
and two-component fermions (/) at finite temperature 
as a function of the boson-fermion (b-f) interaction. In 
weakly coupled b-f mixtures, an induced interaction be- 
tween bosons arising from density fluctuations of fermons 
modifies the critical temperature of the Bose condensate. 
In addition, an induced interaction between fermions due 
to density fluctuations of bosons may lead to superfluid- 
ity of the fermions fs']. On the other hand, in strongly 



coupled mixtures, boson-fermion molecules are formed 
, which may become superfluid . By analyzing the 
realization of internal symmetries in each phase, we clas- 
sify the types of phase boundaries. We then discuss the 
detailed connection with the hadronization phase transi- 
tion in QCD matter at low temperatures. 

We consider a (nonrelativistic) boson-fermion mixture 
with Hamiltonian density. 
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where (p is the boson and ip the fermion field. We label 
the two internal states of the fermions by spin indices a — 
{1, 1}. We focus for simplicity on an equally populated 
mixture of n bosons and n fermions with ni = — n/2. 

The bare boson-fermion coupling g^^ is related to the 
renormalized coupling g^^ and to the s wave scattering 
length by 
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where ei(fc) = k"^ /2mi {i — b, f) is the single-particle ki- 
netic energy, is the boson-fermion reduced mass, and 
A is a high momentum cutoff. We define tq = (2A/7r)~^ 
as a typical atomic scale. We assume an attractive bare 
b-f interaction {g^j < 0), tunable in magnitude, with A 
fixed so that the scattering length Oj^^ can change sign: 



namely. 



g^j,m^/27r for small negative g^^^, while 



'■bf 



tq for large negative g^j,. We keep the bare 
boson-boson and fermion-fermion interactions fixed and 
repulsive {g^^ > 0, gj^>0); thus, the renormalized cou- 
plings g^^ = 47rabb/mj, and g^j, = ATraff/rrij: are always 
positive. These repulsions are required for the stabil- 
ity of the bosons at weak bare b-f coupling and also for 
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the dominance of two-body molecules (&/) over three- 
body molecules (bbf) and (bff) at strong bare b-f at- 
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|: 3nr/2+ (g,, + UbftW/2 < 0. For mjmj 



2 and 



1/3, 



= 10 , the critical values for the instability are 
Boson-fermion mixtures with an attrac- ^f,-BEC ~ 21 at T = and ^ 2.2 near Tc(6-BEC). 

At high T, thermal pressure stabilizes the system. 

A weak b-f interaction also leads to a boson-induced 
attraction between fermions similar to the phonon- 
induced attraction between electrons in metals li. The 
characteristic length of the induced /-/ interaction in fully 
condensed bosons, ay^j, averaged over the Fermi sur- 
face, is a^-^y = (7r/3)^/^(r7itTO^/m^) ln(l + x'^){kp'i]'^)~^ , 
with X — y/Sn/ AkpU^fj. For m^/mj = 2, kpttfhf — 



traction 

five b-f interaction may be realized in three-component 
cold atomic experiments using, e.g., the hyperfine state 
/ = 1, TO/ = 1} of '^^Rb, mixed with the hyperfine states 
|9/2,-5/2) and |9/2, -9/2} of 

We first lay out the phase structure at weak bare 
b-f coupling, where the dimensionless parameter = 



is ^ 1. In the absence of b-f attraction {rj = 



-00) with weak b-b repulsion {n 



1/3, 
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<C 1), boson con- 



densation (&-BEC) occurs below the critical temperature, 
T°{b-BEC) ~ (1 + 1.32ni/3aftj, + • ■ ■)Tq, with the ideal 
BEC transition temperature Tq = 27r[n/C(3/2)]2/3/m^ 
[1^ . A weak b-f interaction induces an attraction be- 
tween bosons via fermion density fluctuations (Fig. [T]) 
given by t/fa/6(p,w;T) = -g2^n(p, c^; T), where n(p, w; T) 
is the one-loop fermionic polarization at temperature 
T. In the static limit, relevant for elastic scattering of 
bosons. 



n(p,0;r) = -2 



d^q "/(|P + q|) - rifiq) 
(2^)3 efi\p + q\)-efiq)' 
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with n^{p) = i/[e[^/(p)-^/l/T + 1] the Fermi occupa- 
tion. Then, the characteristic length of the fermion- 
induced b-b interaction near Tc{b-BEC) is ajj^fj = 

{mj4:n)Ubfbiy^'2mf,To,0;To) < 0, where we set T = Tq 
and take a typical momentum transfer p equal to the 
thermal momentum ■y/2TOj,To. For m^^/m^ — 2, corre- 
sponding to a ^^Rb-^'^K mixture, k^a^^^ 2± —11 77^^ with 
fcy — (37r^ri)i/'^. Since the induced attraction tames the 
bare b-b repulsion, the transition temperature for 6-BEC 
decreases to 



Tc(6-BEC) = [1 + 1.3271^/3(0 
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The replacement a 

der in the b-f interaction; more generally such a replace- 
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a^fb is exact to second or- 



ment is exact in a large M extension (a^ 



where only bubble summations are important [14I . I15l |. 
Calculation of higher order corrections to rc(5-BEC) 
is beyond our present scope. The Bose gas becomes 
mechanically unstable when the net b-b attraction be- 
comes large enough to overcome the thermal pressure 




FIG. 1: Fermion-induced interaction between bosons to sec- 
ond order in the boson-fermion interaction. The single lines 
denote free fermions and the double lines free bosons. 



-4.61n(l -I- x'^)rj ^. A net attraction, 
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< 0, 



causes fermionic BCS superfluidity (/-BCS) below a crit- 
ical temperature [l3| 
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where — fc^/2my and 7 is Euler's constant. For 
"^fc/"^/ — 2, /-BCS emerges when 77 < 77^"^^^ — 
2.1[ln{l-\-x^)/kpaj:jY^'^. As the b-f attraction increases, 
Tc(/-BCS) increases due to the increase of the net /-/ 
attraction. 

Let us turn to the regime of strong bare b-f coupling 
where 77 is large and negative. Here bound molecules or 
composite fermions, N = (bf), are formed with a kinetic 
mass rrij^ = rrif^ -\- rrij . We estimate the s wave scatter- 
ing length of two TV's of opposite spins from the diagram 
sketched in Fig. [2] Integration over the three interme- 
diate states for this process yields the corresponding T- 
matrix element with zero initial and final momenta, 



TwAr(0,0) = 4 / 

J\c^\< 
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where ejv is the binding energy of the composite fermion. 
The prefactor 4 arises from four possible time orderings of 
the dissociation and recombination, (12;34) as in Fig. [21 
(21;34), (12;43), and (21;43). Using the Schrodinger 
equation, [— V^/2to^ -|- V{r)]iljN{r) = eNfpN{^), where 
V(r) is the bare b-f potential, 7/'Ar(r) the internal wave 




FIG. 2: T-matrix diagram for scattering of composite 
fermions of opposite spins, arising from exchange of the con- 
stituent bosons. The single, double, and triple lines denote 
free fermions, bosons and composite fermions, respectively, 
and the three dashed lines indicate the intermediate states. 
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function of the composite fermion, and ipN (q) its Fourier 
transform, we write the matrix element Ai for dissocia- 
tion of A'^ into b and / as = (e^v — 9^/2TOH)V'Ar(q); 
thus, 

T^^.(0,0) = 2 / ^\i,^[q)\\e^-q^/2m,). (7) 

For a bound state in a zero-range potential, 'iI'n{'<^) = 
e ~ /ry/^Tro^ and ejv = — l/2m^a^^, and thus M — 

— y^27r/m^a^ 18]. Equation ([7]) yields the s wave scat- 
tering length for composite fermions of opposite spins for 
negative and large rj 

««« = lfWo,o)^-^a,,r. (8) 

Here 0.76 < r(a^y/ro) < 1 in the interval 1 < a^^/ro < 
oo. This result is the same in magnitude but oppo- 
site in sign from the scattering length between difermion 
molecules in the same approximation 10]. For m^/m^ — 
2, we obtain fc^a^^ ~ — 7.0r/]7y]. The above estimate 
of the A^-A^ scattering length is the leading order term 
in the larg e Af extension of the present boson-fermion 
model [l5|. For finite A/", it is corrected by (i) other 
recombination diagrams [19], (ii) short range b-b and /- 
/ repulsions, and (iii) possible contributions of trimers 
such as (bff) and (bbf) [Ilj. We leave the calculation 
of such corrections for the future and treat Eq. ([8]) as 
a first estimate. Equation ([8]) implies that the low en- 
ergy effective interaction between composite fermions in 
the spin-singlet channel is weakly attractive; the stronger 
the bare b-f coupling, the weaker the A^-A^ interaction. 
Such an effective attraction causes composite fermions to 
become BCS-paired (A^-BCS) below a transition temper- 
ature, 

7 / 9 \ ^/'^ 

Tc{N-BCS) = —(-] e„e'^/2fe^''«iv . (g) 

where e^^ — k'^^/2m^ is the Fermi energy of the A^. 

The above analyses suggest possible phase structures of 
boson-fermion mixtures in the T-rj plane, shown in Fig.[3l 
The system at low temperature for rj = oo is a, mixture 
of the 6-BEC and degenerate unpaired fermions. As 77 
decreases from the right, the size of the 6-BEC region 
shrinks, Eq. @. The /-BCS state emerges at T = T^if- 
BCS), Eq. (O. Both effects are caused by the induced 
attractions. On the other hand, for strong bare b-f cou- 
pling (77 large and negative), the system at low tempera- 
ture is in the A'^-BCS state. As increases from the left, 
the size of A^-BCS region increases according to Eq. 

At intermediate bare b-f coupling (77 ^ 0) where a 
transition from the 6-BEC (and at low T with coexisting 
/-BCS) phase to A-BCS takes place, the phase diagram 
has complex structures depending on the relative magni- 
tudes of ^j^, Qjj^ and g^j. When the intrinsic 6-6 repul- 
sion is weak, uniform 6-BEC at low temperature becomes 
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FIG. 3: Schematic phase structures of boson-fermion mixtures 
with a tunable boson-fermion interaction, (a) Weak g^^; in 
the regime labeled "collapsed" the effective boson-boson in- 
teraction is negative, (b) Strong g^^; a possibly "collapsed" 
region is not shown. 

unstable at 77 = t^^-bec ^ q ^^^^ ^^-^ attraction, 

while at high temperature the system becomes a stable 
normal gas of bosons and fermions because of thermal 
pressure [l^. This situation is indicated by the "col- 
lapsed" region at 77 ~ in Fig. ^a). For strong bare 
6-6 repulsion, the weak coupling formula (|3|) is no longer 
valid, and the 6-BEC phase may possibly survive into the 
region 77 < until the bosons become bound in composite 
fermions as shown in Fig. 3(b). To determine the physi- 
cal structure of the "collapsed" region requires evaluating 
the free energy to higher order. 

To identify the precise phase boundaries in the region 
7; ~ 0, we need to solve the system beyond the weak or 
strong coupling regimes studied here. Rather, we classify 
the phases of the boson-fermion mixture by their realiza- 
tions of the internal symmetry of the system. We focus 
only on the continuous symmetries here. The Hamil- 
tonian density, Eq. ([T|), has U{l)b x U{l)f^ x U{l)f^ 
symmetry corresponding to independent phase rotations 
of 0, Vt> and V^x. On the other hand, 6-BEC, /-BCS 
and A-BCS break U{l)b, U{l)f^+f^ and U{l)b+^f^+f^) 
symmetries, respectively. Here U{1)a±b denotes an in- 
phase rotation of A and B for (-f), and an opposite- 
phase rotation for (— ). Therefore, in the coexisting 6- 
BEC and /-BCS phase, the symmetry breaking pattern is 
U{\)bxu\l)f^xU{l)f^ -^U{l)f^-f^. On the other hand, 
the A^-BCS phase has the symmetry breaking pattern, 
U{l)b X U{l)f, X U{l)f, ^ Uil)t-[f,+f,) X C/(l)/,-/,. 
The difference of unbroken symmetries between these 
two phases implies the existence of a well-defined phase 
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boundary, as indicated in Fig. [SJd, in contrast to the con- 
tinuous crossover from BEC to BCS in a two-component 
Fermi system (20j . 

Interesting problems remaining for further research on 
the phase structure of the mixtures include understand- 
ing the A^-A^scattering length beyond leading order in 
large A/", the quantitative description of the phase struc- 
ture in the intermediate coupling regime, and extensions 
to spin-dependent b-f interactions and unequal popula- 
tion of bosons and fermions. 

The phase structures we find for boson-fermion mix- 
tures of cold atoms display features of those in two-flavor 
QCD with equal numbers of up (u) and down (d) quarks 
with three colors (R,G,B). The ground state of this sys- 
tem at high density is a 2-flavor color superconductivity 
(2SC) with s wave spin-singlet pairing, e.g., between uR 
and dG, in the color antisymmetric and flavor antisym- 
metric channel. The order parameter for color-symmetry 
breaking, SU{3)c — > SU{2)c, is the diquark condensate 
(63) with the diquark operator 6-y = Uj^ajl-ylaClbffp] 
here i,j are flavor and a,/3 color indices, and C denotes 
charge conjugation. The gap is of order a few tens of 
MeV; the remaining quarks, uB and dB, are unpaired 



and form degenerate Fermi seas [21| . On the other hand. 



the ground state of two-flavor QCD with equal num- 
bers of u and d quarks at low density is nuclear matter 
with equal numbers of neutrons and protons, a super- 
fluid state with a pairing gap of a few MeV (see, e.g., 
[221]): the order parameter for the spontaneous breaking 
of baryon- number symmetry U{\)b is the six-quark con- 
densate {N]Nl) = {{b^ct^,^){h0l^)). 

If we model the nucleon, of radius ^ 0.86 fm, as a 
bound molecule of a diquark (of radius r^-, 0.5 fm) and 
an unpaired quark, we can make the following correspon- 
dence between cold atoms and QCD: h ^ 2SC-diquarks 
and / <-!■ unpaired-quark, N ^ nucleon, b-f attraction 
^ gluonic attraction, 6-BEC ^ 2SC, and A^-BCS nu- 
cleon superfluidity [23]. In particular, the boundary be- 
tween the &-BEC and A^-BCS phases shown in Fig. 3(b) 
is strongly suggestive of that deduced between the color 
superconducting quark phase and the superfluid hadronic 
phase in Refs. [2, IfiJ . 
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